Unzipping Vortices in Type-II Superconductors 
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The unzipping of vortex lines using magnetic-force microscopy from extended defects is studied 
theoretically. We study both the unzipping isolated vortex from common defects, such as columnar 
pins and twin-planes, and the unzipping of a vortex from a plane in the presence of other vortices. 
We show, using analytic and numerical methods, that the universal properties of the unzipping 
transition of a single vortex depend only on the dimensionality of the defect in the presence and 
absence of disorder. For the unzipping of a vortex from a plane populated with many vortices 
is shown to be very sensitive to the properties of the vortices in the two-dimensional plane. In 
particular such unzipping experiments can be used to measure the "Luttinger liquid parameter" of 
the vortices in the plane. In addition we suggest a method for measuring the line tension of the 
vortex directly using the experiments. 
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I. INTRODUCTION 

The competition between thermal fluctuations, pin- 
ning and interactions between vortices leads to many 
novel physical phenomena in type-II high-temperature 
superconductors!. Examples include the melting of the 
Abrikosov flux-lattice into an entangled vortex-liquid- 
and the proposed existence of low temperature Bose- 
glass"^, vortex glass- and Bragg glass^ phases. 

Many experimental probes have been used to study 
these phenomena. They include decoration, transport 
and magnetization measurements, neutron scattering, 
electron microscopy, electron holography and Hall probe 
microscopes. More recently it has become possible to ma- 
nipulate single vortices, for example using magnetic force 
microscopy (MFM)^. These can, in principle, measure 
directly many microscopic properties which have been 
up to now under debate or assumed. The possibility of 
performing such experiments is similar in spirit to single 
molecule experiments on motor proteins, DNA, and RNA 
which have opened a window on phenomena inaccessible 
via traditional bulk biochemistry experiments^. 

In this spirit Olson-Rcichhardt and Hastings— have 
proposed using MFM to wind two vortices around each 
other. Such an experiment allows direct probing of the 
energetic barrier for two vortices to cut through each 
other. A high barrier for flux lines crossing has important 
consequences for the dynamics of the entangled vortex 
phase. 

In this paper we introduce and study several experi- 
ments in which a single vortex is depinned from extended 
defects using, for example, MFM. A brief account of the 
results can be found in Ref . [§| . First we consider a setup 
where MFM is used to pull an isolated vortex bound to 
common extended defects such as a columnar pin, screw 
dislocation, or a twin plane in the presence of point disor- 
der. Using a scaling argument, supported by numerical 
and rigorous analytical results, we derive the displace- 
ment of the vortex as a function of the force exerted by 
the tip of a magnetic force microscope. We focus on the 



behavior near the depinning transition and consider an 
arbitrary dimension d. We argue that the transition can 
be characterized by a universal critical exponent, which 
depends only on the dimensionality of the defect. We 
show that unzipping experiments from a twin plane di- 
rectly measures the free-energy fluctuations of a vortex 
in the presence of point disorder in d = 1-1-1 dimensions. 
To the best of our knowledge, there is only one, indirect, 
measurement of this important quantity in Ref. [IQ] . The 
form of the phase diagram in the force temperature plane 
is also analyzed in different dimensions. Related results 
apply when a tilted magnetic field is used to tear away 
vortex lines in the presence of point disorder, which was 
not considered in earlier work on clean systems ^^i. Fur- 
thermore, we show that a possible experimental applica- 
tion of the scaling argument is a direct measurement of 
the vortex line tension in an unzipping experiment. As 
we will show in this paper, in a system of finite size, the 
displacement of the flux line at the transition depends 
only on the critical force exerted on the flux line by the 
MFM tip, the flux line tension and the sample thickness. 
Thus unzipping experiments can provide valuable infor- 
mation on the microscopic properties of flux lines. 



Next we consider a setup where a single vortex is pulled 
out of a plane with many vortices. It is known that the 
large-scale behavior of vortices in a plane is characterized 
by a single dimensionless number, often referred to as the 
Luttinger liquid parameter due to an analogy with bosons 
ind = 1 + 1 dimensions. Wc show that experiments which 
unzip a single vortex out of the plane can be used to 
directly probe the Luttinger liquid parameter. We also 
discuss the effects of disorder both within the defect and 
in the bulk with the same setup. 
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II. UNZIPPING A VORTEX FROM A DEFECT 
A. Review of clean case 

We begin by considering the unzipping of a vortex from 
an extended defect in a clean sample. For a columnar 
defect the system is depicted in Fig. [TJ At the top of the 
sample the MFM applies a constant force f which pulls 
the vortex away from the defect. We assume that at the 
bottom of the sample the vortex is always bound to the 
defect at a specific location. This assumption will not 
influence the results since below the unzipping transition 
the flux line is unaffected by the boundary conditions at 
the far end of the sample. 




FIG. 1: A MFM tip applies a constant force f which pulls the 
vortex away from the defect. The configuration of the vortex 
is represented by r(r). We assume throughout that the vortex 
is always bound to the defect at the bottom of the sample so 
that r(r = 0) = 0. 

In the absence of external force, and for an external 
field aligned with the defect, the appropriate energy for 
a given configuration r(r) of the vortex is given by^: 



L 

dT 



7 



(a.r(r))2 + ^(r(r)) 



(1) 



Here 7 is the line tension and L is the length of the 
sample along the r direction. The vector y{t) represents 
the configuration of the vortex in the d dimensional space 
and V{y) is a short-ranged attractive potential describing 
the d'-dimensional extended defect (in Fig. [T] c? = 3 and 
d' ~ I). The effect of the the external force, exerted 
by the MFM, can be incorporated by adding to the free 
energy the contribution 



Fi = -f • r(L) 



f • (?T-r(r) dr 



(2) 



where we have used t{t = 0) = 0. Here f stands for 
the local force exerted by the MFM in the transverse 



direction. The free energy of a given configuration of the 
vortex is given by 



F(r) = Fo(r) + Fi(r) 



(3) 



The problem, as stated, has been studied first in the con- 
text of vortices in the presence of a tilted magnetic field^^ 
and the results have been applied to the related problem 
of DNA unzipping^-. 

We note that a similar setup can be achieved by using 
a transverse magnetic field instead of the external force. 
See Fig. Indeed in the free energy ^ the exter- 




FIG. 2: Same as in Fig.[T]but with a transverse magnetic field 
instead of the MFM force tearing the flux line away from a 
defect. 

nal force couples to the slope of the flux line drV in the 
same way as the external magnetic field does^. The only 
difference between the two setups is that there are now 
equal and opposite forces acting on the top and bottom 
ends of the sample. However this difference is important 
only in short samples, where the two ends of the flux line 
are not independent from each other. 

In this paper we focus on the thermal average of dis- 
tance of the tip of the vortex from the extended defect 
(xm(r = L)). This quantity is related to the thermal 
average of the length of the vortex that is unzipped from 
the defect, (rm), through {xm) = f{'''m)ll- Here and 
throughout the paper (. . .) denotes a thermal average 
while an overbar denotes an average over realizations of 
the disorder. 

As stated above the universal behavior of (rm) (or 
equivalently (xm)) within this disorder-free model have 
been derived previously. Here we sketch two approaches 
which will be generalized to samples with quenched dis- 
order in the rest of the paper. 

In the first approach, instead of directly summing over 
configurations of the vortex we perform the sum in two 
parts by dividing the vortex into bound and unbound seg- 
ments. The unbound segment starts at the point where 
the vortex departs the defect without ever returning to 
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hit it again up to the top of the sample. Using Eq. ^ 
it is straightforward to integrate over vortex configura- 
tions to obtain for the partition function of the unzipped 
segment 



(4) 



so that the free energy associated with this conditional 
partition function is 

Tu{Tm) = In Zu{Tra) = -/V„/27 , (5) 

where (3 is the inverse temperature. Henceforth in this 
paper we set (3=1, which can be always achieved by 
appropriate rescaling of the energy units. Even though 
the above sum also runs over configurations which re- 
turn to the defect it is easy to verify that these configu- 
rations give rise to exponentially small correction in the 
Tm- Equation ([5]) implies that as the force, f, increases 
the free energy density of the unzipped portion of the 
vortex decreases. In contrast, the free energy density of 
the bound part is, clearly, independent of the force and 
given by Th(Tm) = Vq{L — Tm), where Vq is the free en- 
ergy per unit length of a bound vortex and L is the length 
of the sample along the defect. The vortex will be un- 
zipped when f = fc ^ ^y2^\Vo\ such that the free-energy 
densities of the bound and unzipped states are equal. 

In this representation the total free energy of the vortex 
is given by 



(6) 



The unconstrained partition function of the model is 
given by 



Z = 



dT„ 



-(/'-r)r™/27 



(7) 



Since both results are independent of the dimensionality 
of the defect (columnar or planar) near the transition one 
always finds in the i — )■ oo limit 



1 



(8) 



with ly = 1. Note, that it can easily be seen that ap- 
proaching the transition from above the average length 
of the vortex which is bound to the defect, {{L — Tm)), 
diverges in the same manner^^. 

An alternative approach, which will also be useful 
in this paper, uses the mapping of the problem to the 
physics of a fictitious quantum particle^'^. The contribu- 
tion of the external field f to the free energy now man- 
ifests itself as an imaginary vector potential acting on 
the particle in d — 1 dimensions (with the r axis acting 
as a time direction). Explicitly, usin g th e standard con- 
version from path-integrals (see Ref. for details) one 



finds that the problem can be described in terms of a 
non-Hcrmitian Hamiltonian: 



n 
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-f-p + V{r) , 



(9) 



where p = iV is the momentum operator. In this lan- 
guage the vortex is bound to the defect as long as there is 
a bound state in the Hamiltonian. As mentioned above 
ii is equivalent to a constant imaginary vector poten- 
tial. This analogy makes it apparent that solutions of 
the non-Hermitian problem can be related to those of 
the Hermitian Hamiltonian (where one sets f = 0) by 
an imaginary gauge-transformation^. In particular the 
left '0^(r, f) and the right ip^{r,i) eigenf unctions of the 
non-Hcrmitian problem can be obtained from those of 
the Hermitian problem, -0„(r, f = 0), using 

V'^(r,f) = ZiVn(r,f = 0) 

V'^(r,f) = Vn(r,f = 0)^^-1 , (10) 



where 



u- 



-fr 



(11) 



The universal behavior of Tm at the transition, Eq. 
was obtained in Ref. [ll| by noting that 



JxTp^{r)dv 
/V^(r)dr 



1 



/c-/' 



(12) 



where ipni^) °^ at long r. We note that the imag- 

inary gauge transformation is justified only at / < fc- 
Otherwise the integrals in Eq. p2| formally diverge. In 
principle this divergence can be fixed by proper treatment 
of the boundary conditions^A. 

Finally we discuss the phase diagram of the model. 
It is natural that as the temperature rises the value of 
the critical force needed for the unzipping decreases (at 
very low temperatures a possible reentrance has been dis- 
cussed in literaturei^ii^) . Here we focus on the existence 
of a critical force for any strength of the attractive po- 
tential. The question is completely equivalent to that 
of the existence of a localized vortex line in the absence 
of the external force. Using the analogy to the quan- 
tum mechanical problem it is well known that in two 
dimensions (corresponding to a three dimensional sam- 
ple) and below, as long as drV{r) < there exists 
a bound state. Therefore, in real samples, as long as the 
potential satisfies this condition there is always a mini- 
mum nonzero critical force required to unbind the flux 
line from the defect. 



B. Scaling arguments for the disordered case. 

1. Critical Behavior 

We now consider the effect of point disorder. In this 
case the free energy of a given vortex configuration with- 
out the contribution from the force exerted by the MFM 
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is given by 



L 

dT 



1 



{^rv{T)f + V{v{T))+^,{v{T),r) . (13) 



The frozen point disorder fi(r,T) is assumed to be un- 
correlated and Gaussian distributed with JI = and 
/i(r, T)/i(r', r') — a6{r — r')S{T — t'). As mentioned above 
the overbar denotes an average over reahzations of the 
disorder and the (...) denotes averaging over thermal 
fluctuations. The contribution from the force exerted by 
the MFM retains the same form as in Eq. 

The free energy, Eq. (fT^ . without the contribution 
from the external force has been studied previously with 
and without the presence of an extended defect. With- 
out the defect it is well known that if one fixes one end of 
the vortex at the bottom of the sample the mean-square 
displacement, {x'^{t)) of the vortex after traveling a dis- 
tance T into the sample behaves as (x^(t)) — Bt'^^'^'^\ 
where C,{d = 2) = 2/3 and (^{d — i) k, 0.61 is the wan- 
dering exponent^^. Moreover, for a given realization of 
disorder there is typically a single dominant trajectory of 
the vortex and realizations with two competing minima 
are very rare. Finally, it is known that the free energy 
fluctuations scale as 



5Ti 



(14) 



with uj{d = 2) = 1/3, u}{d = 3) « 0.22. The exponents 
uj and C are not independent. They satisfy the scaling 
relation: oj = 2^—1. Note that as long as there is no 
defect the free energy fluctuations are unaltered even in 
the presence of a force acting on the vortex because this 
force can be simply gauged away. 

The universal properties of the unzipping transition 
with point disorder can be analyzed using a simple scaling 
argument adapted from for the unzipping of DNA 
(for an abbreviated account, see Ref. [§]). 

Following the discussion in the previous section, which 
analyzed the unzipping problem with no disorder, we sum 
over configurations of the unzipped part of the vortex and 
the zipped part of the vortex separately. As before we 
denote the free energy of the unzipped segment by Tu 
and the free energy of the bound segment by Tb- The 
total free energy is then given by 



^{Tm) = Pu{Tm) + J^h{i^ — Tm) ■ 

For large L and t„i we can rewrite 

Tb{L - T,n) ^ Tb{L) - TbiTrn)- 



(15) 



(16) 



For a one dimensional defect, like a columnar pin, J-b{L) 
is a constant equal to the free energy of a flux line com- 
pletely localized on the pin. For higher dimensional de- 
fects like a twin plane, there is a subtlety. To see this, 
consider a zero temperature case first. Then the parti- 
tion function is determined by the optimal trajectory of 
a flux line (in the presence of point disorder) starting at 



T = and ending at t = L — Tm- However, this trajec- 
tory might not be a part of the optimal trajectory going 
all the way from r = to r = L. In this case clearly 
the relation docs not hold. Since this relation does 
not hold at zero temperature, it should not hold at finite 
temperatures as well. However, in Ref. [18] it was ar- 
gued that such situations (where the optimal trajectories 
do not coincide) are very rare and can be neglectecti^. 
Thus ignoring the constant term J-b{L) we can rewrite 
Eq. HSl) as 



(17) 



We can identify three contributions to the free energy 
above. The first is due to the average free energy dif- 
ference between a vortex on the defect and in the bulk 
of the sample. Close to the transition, similarly to the 
analysis in the clean case, it is linear in Tm and behaves 
as a(/c — /)t„i, where a ~ fc/l is a positive constant. 
The second contribution 5Tb, comes from the free-energy 
fiuctuations arising from that part of the point disorder 
which is localized on or near the defect. For a colum- 
nar defect of dimensionality d' = 1 this contribution is 
due to the sum of the fluctuations in the free energy 
about the mean. The central limit theorem implies that 

1 /2 

it behaves as Tm at large Tm- For d' > 1 this result is 
modified and one can use known results for the free en- 
ergy of a direct path in a random media (see Eq. (jl4p '). 

which leads to 5Tb , where d' is the dimension- 

ality of the defect^. Finally, there is a contribution to 
the free energy fluctuations from the interaction of the 
unzipped part of the vortex with the bulk point disorder, 
5Tu- This contribution behaves similarly to 5Tb with a 
different bulk exponent: 5Tu oc Tm '^\ where d> d' \s the 
dimensionality of the sample. Collecting all three terms 
gives: 

T{Tm) = a(/c - f)Tm + 5Tb{Tm) + 5Tu{Tm) ■ (18) 

As discussed above, u}{d) has been studied extensively 
in the past and it is well known that w{d') > w{d) for 
any d' < S^. Therefore, disorder on or close to the de- 
fect controls the unbinding transition for any dimension 
when Tm is large and the problem is equivalent to unzip- 
ping from a sample with point disorder localized on the 
defect. In particular, this result implies that unzipping 
from a columnar defect with point disorder in the bulk 
is in the same universality class as unzipping of a DNA 
molecule with a random sequence of base pairing ener- 
gies. In fact, point disorder is likely to be concentrated 
within real twin planes and near columnar damage tracks 
created by heavy ion irradiation and near screw disloca- 
tions, strengthening even more the conclusions of this 
simple argument. 

The disorder averaged partition function is dominated 
by the minimum of the free energy and thus by configu- 
rations with STb{Tm) < 0. Using T{Tm) ~ a{fc - f)Tm - 

bm^'^'^ \ where 6 is a positive constant, the partition func- 
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tion 



can be evaluated using a saddle-point approximation. We 
then find the value of Tm at the saddle point satisfies 
a(/c - /) = w(d')6T™^'*'^"\ Therefore 



1 



(/c - /) 



with 



Here /c is the critical force and (xm) is the displacement 
of the MFM tip at the transition. To derive Eq. ([23]) 

(19) \)S,&^ the fact that for the unbound segment (xm) = 

We emphasize that the unzipping transition is first or- 
der both in the clean and disordered cases. Indeed, the 
unzipping occurs only at the boundary and does not af- 
fect total free energy of the flux line in the thermody- 
namic limit. Nevertheless, similar to wetting phenomena 

(20) near first order transitions, this transition possesses scal- 
ing properties characteristic of second order transitions 
like diverging correlation lengths, finite size scaling etc. 



1/ = 2 for d' = 1 
V = 3/2 for d = 2. 



(21) 



The result for the columnar defect (d! — 1) agrees with 
known results from DNA unzipping—. 

To check the scaling argument we have preformed nu- 
merical simulations in d = 1+1 dimensions and have been 
able to solve analytically the closely related problem of 
vortex unzipping from a wall in d = 1-1-1 dimensions. 
We have also studied the simplified problems of unzip- 
ping from a d' = 1 and d! = 2 dimensional defect with 
disorder localized on the defect. While the d! — 1 prob- 
lem was solved previouslyi2i22, below (in Sec. Ill C|) we de- 
scribe a new method applying the replica method. Using 
some approximations, this approach can be generalized to 
a. d! = 2 dimensional defect. In the next sections we de- 
scribe all these results which support the simple scaling 
argument presented above. 



2. Finite size scaling and determination of the flux line 
tension. 

Another consequence of the result is the possibility 
to perform a finite size scaling analysis^i. In particular, 
near the unzipping transition the unzipping length (Tm) 
should be a function of the reduced force e cx (/c — /) and 
the system size L. If the scaling ansatz (PO)) is correct 
then we must have 



(22) 



where g{x) is some scaling function. Quite generally we 
expect that when a: 3> 1 finite size effects are unimpor- 
tant and g{x) ^ so that we recover the scaling ((20)) 
and at a; ^ 1 we have g{x) — > go, where go is a con- 
stant. Note that the constant go does not depend on the 
system size. This constant is a universal number of the 
order of one, which depends on the dimensionality of the 
defect. As we find below for the unzipping from a colum- 
nar pin go = 0.5 and for the unzipping from a twin plane 
go ~ 0.7. Relation ([22)1 can be used to extract the line 
tension 7 through: 



7 = fcgo- 



L 



(23) 



3. Phase diagram 

Next, we use scaling arguments to consider the behav- 
ior of the critical force as the strength of the disorder is 
varied. In particular we focus on the existence of a crit- 
ical strength of the disorder beyond which the flux line 
spontaneously unzips without any external force. The 
disorder induced unbinding transition of the vortex from 
the defect in the absence of the force has been consid- 
ered previousl^ r^^'^'^i^^'^^i^^'^'^ . Below we extend a scal- 
ing argument presented first by Hwa and Nattermann in 
Ref. [28| for a columnar pin in arbitrary dimensions to 
include also planar defects. 

Assume that the vortex is localized within a distance 
l± from a columnar pin or a twin plane. Then, it consists 
of uncorrelated segments of length Zii related to l± via 



(24) 



where C is the wandering exponent defined above. Each 
of these segments has a free energy excess of order l'^ 
higher than the energy of the delocalized vortex. The free 
energy cost per length of localization therefore scales as 
Zjj'^^ ^ l^^ . Clearly, a strong enough pinning poten- 
tial gives rise to a constant energy gain per unit length, 
which suppresses the random energy cost of localization 
(note that w < 1 in any dimension). 

So far we established that a localized phase can ex- 
ist. Now let us consider perturbative effects of a weak 
attractive potential and ask whether the vortex imme- 
diately becomes bound to the defect. The free energy 
gained in the presence of the defect, SF, can be inferred 
by perturbations in the strength of the defect pinning 
energy V. Assume that one end of the vortex is held on 
the defect. Then the energy gain due to the attractive 
potential by the defect is associated with the the return 
probability of the flux line back to the defect. Since the 
root mean square displacement behaves as In the return 



probability behaves as 1 1 



-{d-d'K 



The free energy gained 



by hitting the defect SF therefore scales as /|| . 

To determine if the pin is relevant one has to compare 
this energy to the intrinsic variations in the free energy, 
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AF, which scale as I' 



j{d) 



This yields 

n — oc If, 

AF II 



with 



e = l-Cid-d')-uj = 2-{d + 2- d')C- 



(25) 



(26) 



When e < the defect potential is irrelevant, i. e. the 
system gains more energy by minimizing returns to the 
defect, while if £ > it is relevant, i. e. long excm'sions 
are energetically costly. 

As mentioned above in d = 3 numerical simulations 
indicate that C ~ 0.6 which gives for the planar defect 
{d' = 2) e w 1/8 and for the columnar pin {d' = 1) 
£ « ^1/2. Therefore, a weak twin plane is always rele- 
vant and the vortex is always bound to it. However, a 
weak columnar pin is irrelevant and then one expects an 
unbinding transition. In d = 2, where there can be no 
twin plane, C = 2/3 and the columnar defect is found to 
be marginal. As argued in Ref. [28j it is in fact marginally 
relevant. 

To summarize this discussion for columnar defects in 3 
dimensional samples we expect there is a critical strength 
of the bulk disorder beyond which the flux line sponta- 
neously unzips even at zero force. In contrast for a planar 
defect in 3 dimensions and for columnar defects in pla- 
nar 2 dimensions superconductors we expect that for any 
strength of the disorder there is a finite non-zero value of 
the force needed to unzip the vortex. 

Next, we will check the scaling (I^H)) and the anticipated 
localization / delocalization behavior for a number of dif- 
ferent situations using both analytical methods based on 
the replica trick and numerical simulations. 



C. Unzipping from a disordered columnar pin 
without excursions. 

We start our quantitative analysis from the simplest 
situation, where one can get exact analytical results. 
Namely, we consider unzipping from a ID pin with disor- 
der localized only on the pin. Additionally we neglect all 
excursions of the vortex line from the pin except for the 
unzipped region. This problem then becomes identical to 
DNA unzipping. In Ref. the authors analyzed this 
problem using a Fokker-Planck approach and indeed de- 
rived f = 2 near the unzipping transition. Here we show 
how the same problem can be solved using the replica 
trick. The solution was sketched in Ref. [23|. Here we 
review the derivation for completeness and provide addi- 
tional details. 

Ignoring excursions of the bound part of the flux line 
into the bulk gives the free energy a particularly simple 
form. We again write it as a sum over the contribution 
from the bound and unbound segments. The bound seg- 
ment contribution is given by J-b{Tm) — Vo{L — Tm) + 
U{t^), where Vb < is the mean value of the 



attractive potential, L is the length of the columnar de- 
fect which is assumed to be very large, and U(Tm) is a 
random G aussian uncorr elated potential with zero mean 
satisfying J7(T„ijC/(r,„J = A(5(rmi - Tm^)- The contri- 
bution from the unzipped part takes the same form as in 
the clean case (see Eq. ([S])). Collecting the two terms 
gives: 



(27) 



As before we work in the units, where fc^T ~ 1. In the 
equation above the deviation from the unzipping transi- 
tion is measured by e = (/^ — p)/2j, where / is the force 
applied to the end of the flux line and fc = ■\/27| Vb| is 
the critical force. In Eq. (j27p we dropped an unimportant 
constant additive term VqL. 

The statistical properties of the unzipping transition 
can be obtained by considering n replicas of the partition 
function Z{t) = cxp(-J^(r))^: 







= [ dTl.. 


. J dTn exp ( 


Jo 





a=l 



(28) 



where the overbar denotes averaging over point disorder. 
The averaging procedure can be easily done for a positive 
integer n. We eventually wish to take the limit n ^ 0. 
First we order the coordinates Tj, where the j'*'* replica 
unbinds from the pin according to: < ri < ''"2 < • • • < 
Tn. Then for r S [0, zi) there are no replicas bound to 
the columnar pin, for r £ [ti,T2) there is one replica on 
the pin until finally for L > r > t„ all n replicas are 
bound to the pin. Using this observation and explicitly 
averaging over the point disorder in Eq. (j28p we arrive 
at: 



L L 

Z" — n\ J dri . . J dTn exp 

T„_i 



A 



(29) 

where we use the convention Tn+i ~ L. The integral 
above is straightforward to evaluate in the L — > 00 limit 
so that 



^LA/2 J_ TT 1 
X X 1 _ If 



„n^LA/2 { r(l-f l/K„-n) 

A J r(i-fi/K„) 



(30) 



where e„ = e -I- An and k„ = A/2e„. The exponential 
prefactor is an unimportant overall contribution of the 
whole columnar pin while the rest of the expression is 
the [L independent) contribution from the unzipped re- 
gion. Interestingly the restricted partition functions for 
the unbinding problem from a hard wall (with no exter- 
nal force) and for the unzipping from a 1 dimensional pin 
are identical and thus there is equivalence between the 
two problems (see Ref. [29| for more details.) 
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The disorder-averaged free energy is given by the hmit 
y = -hm„^o(^- l)/n ia]. With the help of Eq. ^ 
one obtains 



J^ = ln(eK) + 1'(l/K), 



(31) 



where 5* (a;) is the digamma function and k — /S./2e. The 
unzipping transition occurs at e = or equivalently at 
K ^ oo. The expression (j3ip is identical to the one found 
in Ref. [l^l using a Fokker-Planck equation approach, 
supporting the validity of the analytic continuation in n 
for this particular application of the replica calculation. 
It is easy to see that this free energy yields 



1 



(32) 



where ^'("^(a;) stands for the n-th derivative of the 
digamma function. The expression above predicts a 
crossover from (Tm) ~ 1/e for k <ti I (far from the transi- 
tion) to (Tm) ~ n/e — A/e^ for k ^ 1 (close to the tran- 
sition) similarly to the unzipping from the wall problem 
analyzed above. Also, it is easy to check that 



1 



(33) 



Here there is a crossover from w « for k ^ 1 to 
ui « 2K/e^ = A/e'^ for K ^ 1. As has been noted in 
the context of DNA unzippingi^ \/w/{Tm) changes from 
being of order unity for the weakly disordered k <C 1 case 
to ~ e^/^ for K ^ 1. Thus for k ^ 1, close to the unzip- 
ping transition, thermal fluctuations become negligible 
and one can work in the zero temperature limit. 

The simplicity of the problem also allows finding the 
higher moments of the distribution. Here we evaluate the 
second moment, which gives the width of the distribution 
of (tto) due to different disorder realizations. Note that 
since the order of averaging over thermal fluctuations and 
disorder is important this quantity can not be extracted 
directly from Eq. l|33|) . To proceed we consider the gen- 
erating function, yV„(ej ) defined by 



dri 



The second (and similarly the higher) moments can be 
found by differentiating yV„ with respect to e^: 



lim 



1 1 ^ d^Wniej) 



(34) 



Upon evaluating the integral, we find 

1 



Wn(e,) = n 



(35) 



and correspondingly 



= lim — —, , . . 



(36) 



This double sum can be calculated using a trick similar 
to the one described in Ref. 1311: 



2k' 



dxdy- 



1 



ye" 



x>y>0 



— M,W{1/k) (C + *(1/k)), 



(37) 



where C « 0.577 is Euler's constant. In the limit of weak 
dis order or high temperature k <C 1 , not surprisingly, we 



get (r^ 



2/e , which agrees with the Poissonian statis- 



tics of with an average given by (r™) — 1/e. In the 
opposite limit 1 one finds (t^) = 4K^/e^. Note that 
(Tm) — 't/e, thus the relative width of the distribution 
i^Tm/ (Tm)), defined as the ratio of the variance of the 
unzipping length Tm to its mean is larger by a factor of 
•\/3 than that in the high temperature regime. The distri- 
bution thus becomes superpoissonian at large k. In fact, 
in the limit k — > oo one can derive the full distribution 
function -PK^oo(''ni) using extreme value statistic s ^^i'^^ : 



VK^ooiTm) ~ e/KG{Tm £ / k) 



(38) 



G{x) = e"^/-* - ierfc(VJ/2), (39) 



with 



where erfc(a;) is the complimentary error function. It is 
easy to check that this distribution indeed reproduces 
correct expressions for the mean and the variance. We 
emphasize that while the thermal fluctuations of the un- 
zipping length become negligible near the transition, the 
fluctuations due to different realizations of point disorder 
are enhanced and lead to a wider-than-Poissonian distri- 
bution of Tm- 

To check these results and uncover subtleties that 
might arise in experiments, we performed direct numeri- 
cal simulations of the partition function of the free energy 
(P7)) . For this purpose we considered a discrete version 
of the problem where the partition function is 



Z 



(40) 



Here U{jni) is the random potential uniformly dis- 
tributed in the interval [~Uq,Uq\ so that the disorder 
variance is A = U'{mi) = Uq/3. For the simulations 
we choose e = ln(1.2) - 0.18 « 0.00232 and Uq = 0.3, 
which gives A = 0.03, k « 6.46 an d acc ording to both 
Eq. ((32|) and numerical simulations (rm) ~ 2860. Then 
we computed 6Tm/ (Tm) using both Eq. ([37]) and perform- 
ing numerical simulations. For the chosen parameters the 
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equation ([37)) gives (5t,„/ (r„ 
simulations yield i5rm/ (rm) 



i 1.68, while the numerical 
1.67. Clearly the results 



are very close to each other and the small discrepancy 
can be attributed to the dis creti zation error. In Fig. [3] 
we plot dependence of 5Tm/{Tm) vs. system size. It is 
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FIG. 4: Dependence of the relative width of the distribution 
on the disorder parameter n. The two curves correspond to 
different averaging over temperature and disorder (see text 
for details) . The horizontal line at \/3 denotes the asymptotic 
value of both 5rm and (Sr^ at k — > oo 



FIG. 3: D ependence of the relative width of the distribution 
5Tm/ ^Tm) ou the System size. Symbols correspond to the ac- 
tual data, the solid line is the guide to the eye, and the dashed 
line corresponds to the replica result in the thermodynamic 
limit. 

obvious from the figure that in the thermodynamic limit 
L — > oo the replica result is in excellent agreement with 
numerical simulations. We mention that numerical simu- 
lations of 5Tm show very strong finite si ze effe cts. There- 
fore one has to go to very large L > 50 (rm) in order to 
approach the thermodynamic limit for the width of the 
distribution. 

Depending on the system the quantity (r,-^) is not al- 
ways experimentally accessible. For example, in the un- 
zipping experiments it is easier to measure thermal aver- 
age, (tto), in each experimental run. We note that this 
quantity has sample to sample fluctuations only due to 
the presence of disorder. Then t he vari ance of the dis- 
tribution will be characterized by (rm)^. The difference 
between the two expectation values is given by w found 

in Eq. ([55]) . Defining (5r^)^ = (Tm)^ — (Tm) and us- 
ing Eqs. dSZl) and §^ we find that 6t^/{t„,) « 
in the weak disorder limit [k <C 1) and St^/ (rm) ~ 
•\/3— 1/{^/3k) in the opposite limit k ^ 1. We plot both 
St^ and 6Tm. versus the disorder parameter k in Fig. U) 
The same issue of importance of the order of thermal 
and disorder averaging appears in the calculation of the 
higher moments of r^, becoming irrelevant only in the 
limit (k oo), which effectively corresponds to the zero 
temperature case. 



Before concluding this section let us make a few re- 
marks about the rebinding transition, i.e., the rezipping 
that occurs with decreasing force. One can consider a 
similar setup with a lower end of the flux line fixed at 
the bottom of the columnar pin and the top end is pulled 
away from the pin with a force /. However, now we will 
be interested in / > fc- Then clearly most of the flux 
line will be unzipped from the pin except for a portion 
near the bottom end. If / is very large, the length of the 
bound segment fm near the sample boundary is small. 
However as / decreases and approaches fc from above, 
the length of this segment increases and finally diverges 
at the transition. This rebinding transition can be de- 
scribed in a similar spirit to the unbinding. For example 
instead of the free energy (|27|) one has to deal with 



Hrm) = W,r.+ r dT'^U{T'J. (41) 

Jo 



As we already noted the free energies ([?7|) and (|1T|) are 
equivalent up to an unimportant constant equal to the 
total disorder potential of the pin: dr^J7(r^). We 
conclude that the unbinding and rebinding transitions 
for a single flux line on a disordered columnar pin are 
identical. In other words, statistical properties of Tm for 
a given f = fc — Sf are identical to those of tv„ for / = 
fc + Sf. 
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D. Unzipping from a planar defect without 
excursions. 

We now generalize the ideas of the previous section to 
the more complicated problem of unzipping of a single 
flux line from a disordered twin plane. As before we 
ignore excursions out of the plane for the bound part of 
the flux line. Let us consider the rebinding transition 
flrst. That is we assume that / is slightly greater than 
fc and we study the statistics of the bound part of the 
flux line. We again assume that the flux line is pinned at 
the bottom of the plane (r = 0) and unbinds for r larger 
than some fm- 

The point disorder potential now depends on the two 
coordinates r and z spanning the twin plane. Using Eq. 
([5|) the partition function reads: 



Z = J dfm J -Dz(T')exp 



Vf„ 



(42) 



where y < is the mean attractive potential of the 
twin plane and we have dropped the unimportant L- 
dependent factors. As before, we assume a Gaussian 
random noise with zero mean and 



/l<(ti, zi)/z(r2, Z2) = aS{Ti - T2)S{zi - Z2). 



(43) 



We also introduce e — —p/{2'y) — V . Note that for the 
rebinding transition e < 0. After replicating the partition 
function and averaging over point disorder we flnd 



L L 

n\ I df„ I dfn-i 



dTi / Dzi(T[)...DZnK) 







exp 



n ^ 
^ ef^ + / dT'^Ca[zi{T[), 



(44) 



where we define fn+i = and Ca is the Euclidean La- 
grangian corresponding to the Hamiltonian {TLa) of a 
interacting particles^: 



Ha — 



1 



27 dzj 



13=1 



■ E 

1</3<7<Q 



5{zp-z^). (45) 



Close to the rebinding transition, we anticipate f ^ ^ 00 
and thus the mean separation between the rebinding 
times of different replicas fa and fa-i diverges. There- 
fore the contribution to the partition function coming 
from integration over Tq will be dominated by the ground 
state of configurations with a replicas. In this case we 
can significantly simplify the partition function and eval- 



uate it analytically: 

L L L 

'Z^ = n\ j dfn J dfn-i ■■■ j dfi 

f„ f2 



(46) 



exp 



ETq + {£a ~ £a-l)fa) 



Here £a is the ground state energy of Ha with a sub- 
tracted term linear in a, that just renormalizes fc- Close 
to the transition e is linear in the difference f ^ fc- The 
energy, was computed in Ref. [3T||: 



£a = -^a^ = -^a\ 
Upon integrating over fa one obtains 

1 



n 



n 



^ |e|a - ^a^ ^\ \e\ - ^a^ 

a—l ' ' a— 1 ' ' ^ 



(47) 



(48) 



The product above can be reexpressed in terms of F- 
functions, which in turn allows for a straightforward an- 
alytic continuation to n — > 0: 



r 



(49) 



Using this expression we obtain the free energy and the 
mean length of the localized segment: 



J- = — lim 



Z" - 1 



n^O n 

dT _ 



= lne 



2^- 




, (50) 



(51) 



where as before ^'("''(.t) stands for the nth derivative of 
the digamma function. This expression has the asymp- 
totic behaviors: 



2|e|3/2 



S»|e|. 



(52) 



This scaling confirms the crossover between exponents 
V = \ and V — "^12 for the rebinding transition to a 
two-dimensional disordered plane predicted by the simple 
scaling argument leading to Eq. ([2T|) . 

In a similar way one can also consider an unzipping 
transition with / < /c- One finds an expression for the 
partition function which is identical to ([48]) with the sub- 
stitution ^ — !■ — ^. Note however, that the analytic con- 
tinuation of the product (j49p results in a complex parti- 
tion function and hence a complex free energy. It thus 
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appears that the analytic continuation of the product 
(|48p to noninteger values of n is not unique. One can 
always multiply it by any periodic function of n, which 
is equal to unity when the argument is integer. While we 
were able to find some real-valued analytic continuations 
of Z"^ to negative values of ^, these continuations did not 
lead to physically sensible results. 

Because of the ambiguity of the analytic continuation 
and some approximations used to derive Eqs. (|50p and 
(|5ip we also performed numerical simulations for the vor- 
tex unzipping from a disordered twin plane. 

For numerical simulations we are using the lattice ver- 
sion of the model, where in each step along the r direction 
the vortex can either move to the left or the right one lat- 
tice spacing. Note that because we neglect excursions the 
vortex motion occurs strictly within the plane until the 
vortex is unbound. Then the restricted partition func- 
tion for the bound part of the flux line, Z{x,t), which 
sums over the weights of all path leading to x, t, starting 
at a; = 0, r = satisfies the recursion relational 

Z{x, T + l) = e^(^'^+i) [jZ{x -1,t) + JZ{x + 1, t) 
+ {1-2J)Z{x,t)]. (53) 

We assume that ^{x,t) is uniformly distributed in the 
interval [—Uq,Uq\ implying as before the variance a = 
C/o/3. The variable J controls the line tension. In the 
continuum limit J ^ 1 and Uq <^ I the equation (j53|) 
reduces to the Schrodinger equation: 

dZ , , 

— = -nZ{x,T) (54) 

OT 

with the Hamiltonian given by Eq. ([9|) with j — 2J and 
/ = (there is no force acting on the flux line within 
the plane). We note that even if the parameters of the 
discrete model are not small we still expect that Eq. ([M]) 
remains valid at long length and time scales. However, 
the relation between the microscopic parameters of the 
discrete model and the parameters of the effective coarse- 
grained Hamiltonian ^ is more complicated. 

In our simulations we evaluated numerically the free 
energy of the bound part of the vortex line for each real- 
ization of point disorder and used the analytical expres- 
sion for the free energy of the unbound part, for which 
point disorder can be neglected. The latter is given by 
Eq. ([S])- This free energy is controlled by a single param- 
eter 7^/(27). Use of the analytic result ^ significantly 
simplifies calculations of (t,„) and allows us to perform 
large scale simulations. 

First we verify the scaling (|20p with 1/ = 3/2 at the un- 
zipping transition. To do this we perform standard finite 
size scali ng pr ocedure. In Fig. [5] we show dependence of 
the ratio (t„i) /L on the parameter f'^/{2j) for four differ- 
ent sizes. As we expect from the scaling relation the 
three curves intersect at the same point corresponding to 
the unzipping transition {qq « 0.7). Once we determine 
the crossing point corresponding to the critical force fc 
we can verify the scaling relation with = 3/2. In 



1.0 




FIG. 5: Ratio of the unzipping length (xm) to the system size 
L as a function of f'^ /2"f for different system sizes. Here / 
is the external force and 7 is the line tension of the vortex 
(see Eqs. (|4]and (Hjl)). According to the scaling relation ([22} 
the crossing point corresponds to the unzipping transition. In 
simulations the parameters of the microscopic model (I53p are 
chosen to be J = 0.2, Uo = 2 

Fig. [5] we plot (Tm) / L versus the scaling parameter eL^^'^ 
(see Eq. (HU) with = 3/2 for two different system sizes. 
Clearly the data collapse is nearly perfect, which proves 
the validity of the scaling (|63p with 1/ = 3/2 for the unzip- 
ping of a flux line from a twin plane. The inset shows the 
derivative of (rm) with respect to e. Clearly this deriva- 
tive is asymmetric with respect to e ^ — e, implying 
that there is no symmetry between the unbinding and 
rebinding transitions. This is contrary to the unzipping 
from a columnar pin with no excursions, where such a 
symmetry does exist. 

Next we turn to verifying the analytic prediction for 
(t„i), Eq. ([5T|) . As we argued above the parameter ( 
describing the disorder strength can be easily extracted 
from microscopic parameters of the model only in the 
continuum limit Ug 3> 1, J <C 1. Unfortunately, it is 
not possible to do simulations directly in the continuum 
limit ( J <C 1 and Uq <C 1). Indeed as Eq. ([5T|l suggests 
in order to see the scaling exponent u = 3/2 one needs 
to go to length scales much larger than 1/^, where ^ — 
a^J/12 = If J < 1 and especially Uo < 1 

then one has to simulate extremely large system sizes 
where L is larger than 10^ for Uq = 0.1 and J = 0.1. 
Therefore we perform simulations in the regime where 
J and especially Uq are appreciable. We then regard 
^ as a fitting parameter of the model which should be 
equal roughly to t/g J/36. In Fig. [7] we show results of 
numerical simulation for the rezipping length L — (rm) 
on the detuning parameter e for different system sizes. 
The solid black line is the best single-parameter fit to 
the data using the analytic expression (|?T1) . The fitting 



11 



1.0 



0.8 



0.0 



/.=6400 
■/.= 12800 



-J 




■^0.6- 




E 








Vo.4- 


/ " 
/ ^ 




/ 1 A 




/ 

/ I-" 


0.2- 


/ 1 V 



, 1/v 2 
8L 



FIG. 6: Data collapse of {rm) /L as a function of el}^" with 
the exponent i/ = 3/2 for two different system sizes (see 
Eq. (|22p 'l. The parameters of the model are the same as in 
Fig. [5] The inset shows derivative of {r„i) with respect to e 
for L = 12800. Clearly the scaling function is asymmetric 
with respect to e — e. Thus the unbinding and rebinding 
transitions are not equivalent. 
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see that indeed as decreases the fitting parameter ^ 
becomes closer to the continuum expression. 

While we were not able to derive a closed analytic ex- 
pression for (tto) for the unbinding transition, we per- 
formed numerical simulations. As the inset in Fig. [6] 
suggests the transition is highly asymmetric. In fact this 
asymmetry persists in the thermodynamic limit. In Fig. [8] 
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FIG. 8: Comparison of dependences of -L — (rm) for the re- 
binding transition and (rm) for the unbinding transition on 
|e|. We used the parameters of Fig. [5] with L = 51200. The 
finite size effects are negligible on the scale of the graph. Both 
curves interpolate between l/|e| dependence at |e| 2> ^ and 
C/|e|''/'^ at |e| ^ ^. However, the prefactor C for the un- 
binding transition is about three times larger than for the 
rebinding. 



we plot i— (t„i) for the rebinding transition and (t,„) for 
the unbinding versus |e|. Both curves interpolate between 
l/|e| dependence at weak disorder |e| <C ^ and C/|ep/^ 
dependence at strong disorder |e| ^. However, the 
prefactor C in front of l/|ep/^ is larger for the unzipping 
transition. 



E. Unzipping from a hard wall 



FIG. 7: Dependence of the leng th of the bound part of the 
fiux line to the twin plane L — {jm) on t for the rebinding 
transition. Different curves correspond to different system 
sizes. The solid black line is the best single parameter fit 
using Eq. (|51|l with ^ being the fitting parameter. 



parameter ^ found from simulations is ^ w 0.036, while 
a continuum estimate [/q J/36 gives ^ « 0.089, which 
is very reasonable given that this estimate is valid only 
at [/q ^ 1- We also performed similar simulations for 
t/o = 1-5 and got a very good fit with ([5T|l for ^ = 0.018, 
while the continuum estimate gives ^ w 0.028. We thus 



As the next step we consider unzipping from an attrac- 
tive hard wall in c? = 1-1-1 dimensions with point disorder 
in the bulk. Our method is a straightforward generaliza- 
tion of the Bethe ansatz solution found by Kardar in the 
absence of the external force'^^. The system is illustrated 
in Fig. [9] Here the potential experienced by the flux 
line, V(x)^ has a short ranged attractive part and an im- 
penetrable core at X = 0. While the scaling argument 
is unchanged in this case, this problem has the merit of 
being exactly solvable within the replica approach. Since 
most details of the calculation are identical to those pre- 
sented in Ref. [3l|, here we only outline the solution. 
After replicating the free energy Eq. (fT5)) along with the 
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FIG. 9: An illustration of the setup considered in the Bethe 
Ansatz calculation. The flux line is restricted to the half plane 
and a MFM tip is acting on it at the top of the sample. 



contribution from the external field Eq. ^ and averag- 
ing over the disorder, the replicated sum over all path 
weights connecting points (0,0) and (x,r), Z"{x,t), can 
be calculated from 



(55) 



with the initial condition Z"(a;, 0) = 5{x). The replicated 
system describes n attractively interacting bosons with a 
non-Hermitian Hamiltonian 7i given by 



1 



27 



a=l 

- cr ^ 5{xa - Xfj) - ^an . 



(56) 



a<l3 



In Ref. [31| the problem was solved for / = using 
the Bethe Ansatz. The boundary conditions were that 
the ground state wave function should vanish at large 
X should decay as exp(— Ax) for the particle closest to 
the wall. One then finds that for the permutation P of 
particles such that < xpi < xp2 < ... < xpn the wave 
function for / = is 



/=o ~ exp 



KaXpo 



(57) 



Here = A + 2(a — 1)k, k = (77/2. Taking the 
zero replica limit it was found'^^ that for weak disorder 
{(J'y/2 < A) the vortex is bound to the wall while for 
strong disorder (0-7/2 > A) it is unbound. 

The ground state wave function for the non-zero value 
of the force can be obtained by noting that the non- 
Hermitian term acts like an imaginary vector potential. 
In particular, it can be gauged away when the vortices are 



bound to the wall as discussed in Sec. Ill Al (see Eqs. (fTO]) 
and (jlip). This imaginary gauge transformation gives 



n 



(58) 



^-Z = f/^oexp , 
which implies that the solution is 

■^f = exp ( - ^ kaXpa I , (59) 
a— 1 / 



with ka — X + 2{a ~ 1)k — f . The effect of the force 
is simply to shift all the Kq's by a constant. The aver- 
age localization length (which satisfies near the transition 
{Xm) - fc{Tm)h) is then given by 



1 



Yldxj 



(60) 



where Z„ = YVj^i ^^j"^ fi^j)- Note that the nor- 
malization factor Zn in the equation above is formally 
equivalent to the partition function ([29)) for the unzip- 
ping from a columnar pin without excursions if we iden- 
tify A — K — / with e and k with A/2. This equivalence 
implies that (xm) for the unzipping from a hard wall has 
the same statistical properties as (Tm) for the unbinding 
from a columnar pin (for more details see Ref. [29|). In 
particular, the unzipping problem has a crossover from 
(xm) - l/{fc - /) for A - / » « to (x„-,) ^ l/{f, - /)3/2 
in the opposite limit. 

This example confirms another prediction of the simple 
scaling argument: the critical exponents for the unbind- 
ing transition are determined only by the dimensionality 
of the defect even if the disorder is also present in the 
bulk of the system. 



F. Unzipping from a columnar pin with excursions 
into the bulk. 

In this section we consider the setup similar to 
Sec. Ill CI namely unzipping from a columnar defect in 
d = 1 + 1 dimensions, but allowing excursions of the flux 
line to the bulk (see Fig. fTU]) . Unfortunately there is 
no analytic solution available for this problem. There- 
fore we present only numerical results. As in Sec. Ill Dl 
we consider a lattice version of the model where in each 
step along the r direction the vortex can either move to 
the left or the right one lattice spacing. The attractive 
potential was placed at a; = 0. The restricted parti- 
tion function of this model, Z{x^ t), which sums over the 
weights of all path leading to x, t, starting at a; = 0, r = 
satisfies the recursion relation"^^: 

Z(x,T-Hl) = 4,o(e^- 1)^(0, r) 

_^gM(^,r+i) [j^fzix - 1, t) + Je-'fZ{x + 1, t)] (.61) 
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FIG. 10: A setup illustrating unzipping from a columnar pin 
in d = 1 + 1 dimensions with excursions into the bulk. 
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Similarly to Eq. (|53p we assume that ii(x,t) is uniformly 
distributed in the interval [~Uo,Uo] implying the vari- 
ance a = Uq/3. The variable J controls the line ten- 
sion, V is the attractive pinning potential, and / is pro- 
portional to the external force. In the continuum limit 
J <C 1, / <C 1, and Uq 1, equation ((6T|) reduces to the 
Schrodinger equation: 



a7 



= -nZ{x,T) 



(62) 



with the Hamiltonian given by Eq. ([9]) with 7 = 2 J. 

For the simulations we have chosen particular values 
of J = 0.1 and y = 0.1. As before we work in units 
such that ksT = 1. In the results described below the 
partition function was evaluated for each variance of the 
disorder for several systems of finite width w = 2L,j. av- 
eraging over the time-like direction (typically r ~ 10^ 
"time" steps) with the initial condition Z(0, 0) = 1 and 
Z{x,Q) = for x 7^ 0. 

To analyze the numerics we performed a finite size scal- 
ing analysis. In the spirit of Eq. (PO)) . in the vicinity of the 
transition we expect the scaling form (compare Eq. (|22p ): 



(63) 



where $ is some scaling function. Based on the results 
of previous sections we anticipate a smooth interpolation 
between scaling exponents v = 1 and u = 2 with either 
increasing or increasing strength of disorder at fixed 
L^. To perform the finite size scaling we obtain for each 
value of Lx a value for the exponent v from the best 
collapse of the numerical data of two systems sizes 
and Lx/2. In Fig. [Tl] we plot 1/v as a function of the 
system size Lx- As can be seen the data is consistent with 
V saturating aX v — 2 for large systems. The crossover 
to = 2 is much more rapid if the point disorder is 
enhanced near the columnar pin (see the inset in Fig. llip . 
as might be expected for damage tracks created by heavy 
ion radiation. 

Next, we test the behavior of the critical force as the 
disorder strength is increased. According to our discus- 
sion in Sec. IIIB 3[ we anticipate that in the absence of 



FIG. 11: Effective exponent l/u versus Lx for a fixed strength 
of point disorder a = 0.03. The results are consistent with 
the general argument that this exponent should saturate at 
1/ = 2 as — > 00. The inset shows the same exponent vs ctc, 
the variance of additional point disorder placed directly on 
the columnar pin extracted from two system sizes ~ 600 
and Lx = 1200. It appears that — > 2 as CTc increases. 



an external force the flux line is always bound to the pin 
in 1 + 1 dimensions. This is in contrast with the problem 
of unzipping from the wall discussed in the previous sec- 
tion, where there is a critical strength of the disorder, ctc, 
which leads to an unbinding transition for f — 0. Note 
that the existence of a critical value of the disorder is a 
direct consequence (see discussion in Sec. Ill B 3|) of the 
excursions of the vortex from the defect which, as argued 
above, do not modify the critical behavior of the unzip- 
ping transition. The existence of a critical value of the 
disorder is therefore strongly dependent on the dimen- 
sionality of the problem. 



In numerical simulations for each strength of disor- 
der we determine the critical force plotting the ratio 
(Tm) /Lx for two different sizes and using the scal- 
ing relation (l63|) . Note that this ratio does not depend 
on Lx at f — fc (see also the discussion in Sec. Ill D| ). We 
checked that this is indeed the case. Upon repeating this 
procedure for different disorder strengths we obtain the 
dependence fdUa) which is plotted in Fig. [121 The graph 
suggests that there is no unbinding transition at zero tilt 
at any strength of disorder consistent with the scaling 
argument presented in Sec. Ill B 31 and those of Ref. [2^. 
We point out that the strongest disorder shown in the 
graph Uq = 0.9 required samples quite extended in the 
time- like direction, Lr « 10*. 
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FIG. 12: Critical force for unzipping from a columnar defect 
in 1 + 1 dimensions as a function of the disorder strength. 



III. UNZIPPING A LUTTINGER LIQUID 

We now turn to consider the effect of interactions on 
the unzipping of single vortices. To do this we study a 
system where the vortices are preferentially bound to a 
thin two-dimensional slab which is embedded in a three- 
dimensional sample so that the density of vortices in the 
slab is much higher than in the bulk. Experimentally, 
this setup could be achieved using, for example, a twin 
plane in YBCO or by inserting a thin plane with a re- 
duced lower critical field Hd (with, for example, molec- 
ular beam epitaxy) into a bulk superconductor. The sce- 
nario we analyze is one where a MFM is used to pull a 
single vortex out of the two-dimensional slab (see Fig. 
[T5)) . The physics of the vortices confined to two dimen- 
sions is well understood and is analogous to a spinless 
Luttinger liquid of bosons (see, e.g. Ref. Q)- 

As we show below the dependence of the displace- 
ment of the vortex from the two-dimensional slab on the 
force exerted by the MFM depends on the physics of the 
two-dimensional vortex liquid which resides in the slab. 
Specifically, the critical properties of the unbinding tran- 
sition depend on the "Luttinger liquid parameter" which 
controls the large-distance behavior of the vortex liquid. 
The experimental setup can thus be used to probe the 
two-dimensional physics of the vortices in the slab. 



A. Two-dimensional vortex liquids 

The physics of vortices in two dimensions is very well 
understood. The vortices form a one-dimensional array 
located at position Xiir). The density profile of the vor- 



f 



y 



FIG. 13: Possible experimental setup for studying unzipping 
from Luttinger Liquid. A MFM is used to pull a single vortex 
out of a plane where the vortices are confined. The measured 
quantity is the distance of the pulled vortex from the confining 
plane as a function of the force /. 

tices is then given by 

n(x,T)=.^<5[x-x,(r)] , (64) 

j 

where x and r denote transverse and longitudinal coordi- 
nates with respect to the vortices and i is an index label- 
ing the vortices. By changing variable into the phonon 
displacement field Uj through xj (t) ~ a [j + Uj (t)] , 
where a is the mean distance between vortex lines the 
free-energy of a particular configuration can be written 
as: 

J^o = Y / dxdr [ciiid^u)'^ + C44{druf] . (65) 

Here cn and C44 are the compressional and the tilt mod- 
uli respectively. After rescaling the variables x and t 
according to 

x^ xi — ] , T ( — ] , (66) 

VC44/ Veil/ 

the free energy takes the isotropic form 

^0 = ^1 dxdr [(d^u)^ + [druf] (67) 

with A = a^^ciiC44. The partition function is then given 
by the functional integral 

Z = j Du{x,T)e-'^ , (68) 

with S = Sq = To /T. In the limit of large sample sizes 
in the "timelike" direction one can regard Z as the zero 
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temperature partition function of interacting bosonsi^. 
In this language the imaginary time action can be written 
as 



(69) 



Here we set h — 1 and identified the Luttinger-hquid 
parameter, as 



9 



ttT 



(70) 



The Luttinger-hquid parameter controls the long- 
distance properties of the model. For vortices g it is a 
dimensionless combination of the compressional and tilt 
moduli, the density of vortices and temperature. 

Various properties of Luttinger liquids are well under- 
stood. For example, the correlation function for the den- 
sity fluctuations Sn(x, r) — n{x, t) — no, where tiq — 1/a 
is the mean density, obeys 



{Sn{x,T)6n{0,0)) 



cos (27mox) 
(a;2 + r2)9 



(71) 



There is quasi long-range order in the system and the 
envelope of the density correlation function decays as a 
power law with the exponent depending only on g. As 
we show below, g can be probed by unzipping a single 
vortex out of a plane which contains a (l-)-l)-dimensional 
vortex liquid. 

In what follows we also consider the case where there is 
point disorder present in the sample. The behavior will 
be strongly influenced by the behavior of the vortices in 
two dimensions in the presence of disorder. This prob- 
lem has been studied in some detail in the past (see e.g. 
Ref. [1^ and references therein). Here we briefly review 
features which will be important in analyzing the unzip- 
ping problem. The most relevant (in the renormalization 
group sense) contributions to the action from the point 
disorder is 

SpD = 2 / dxdTR{x,T) cos[2t:u{x, t) + /3{x,t)] , (72) 



where positive (negative) R implies a repulsive (attrac- 
tive) potential between the vortices and the quenched 
random disorder. We assume, for simplicity, that (3{x, r) 
is distributed uniformly between and 2tt and R(x, t) 
has a an uncorrelated Gaussian distribution with the 
variance Ag: 



R{xi,Ti)R{x2,T2) = Ao(5(xi - x2)6{ti - T2) , (73) 

where the overbar, as before, represents averaging over 
disorder. 

To analyze the disordered problem, similar to the sin- 
gle vortex case, we use the replica trick. Then the repli- 
cated noninteracting part of the action becomes 



a,f3 



dxdr 



dua dufj dua dup 



dr dr dx dx 



.9 

(74) 



Here Ua{x,T) is the replicated phonon field and k is an 
off-diagonal coupling which is zero in the bare model 
but is generated by the disorder. It plays the role of a 
quenched random "chemical potential" which is coupled 
to the flrst derivative of the phonon field u. The replica 
indices, a and /3 run from 1 to n and at the end of the 
calculation one takes the limit n — > 0. After replication 
the contribution from the point disorder becomes 



Uf3{x,T)] 



S'p£) = —Ao^^ / / dxdr cos 27T [u a {x,t) 

a, 13 

(75) 

The combined action can be treated within the renor- 
malization group using a perturbation series near 5 = 1 
where a phase transition between a vortex liquid and a 
vortex glass occurs-?^-. By continuously eliminating de- 
grees of freedom depending on frequency and momentum 
within the shell A — JA < ^/uj^^+cp < A, one obtains the 
following renormalization group equations^i^ 



^ = 



dg 

dl 
dA 

'dl 
dn 

m 



2(1 - .g)A - 2CA2 
C^A^ 



(76) 
(77) 
(78) 



Here I is the flow parameter K{1) — Ke^^ . C is a non- 
universal constant which depends on the cutoff A. The 
equations are subject to the initial conditions k{1 = 0) = 
and A(/ = 0) = Aq. Note that the Luttinger liquid 
parameter is not renormalized. Analyzing the flow equa- 
tions it has been shown that in the vortex liquid phase 
(5 > 1) the correlation of the density fluctuation behaves 
in the vortex liquid phase as 



((5n(x,T)Jn(0,0)) 



1 



(79) 



where k is a nonuniversal exponent. In the glass phase 
{g < 1) correlations decay faster than a power law, with 

{5n{x, r)(5n(0, 0)) ~ exp (-(1 - gf In^ ^/ x^ + t"^^ . 

(80) 

In what follows we consider a setup in which a two di- 
mensional array of vortices, whose properties have been 
described above, is embedded in a three dimensional bulk 
sample. As shown below when a single vortex is unzipped 
into the bulk in a clean sample the critical properties of 
the unzipping transition yield information on the prop- 
erties of the two dimensional vortex liquid. In particular, 
they provide a direct measure of the Luttinger-liquid pa- 
rameter. In the same setup in a disordered sample we 
will show that the critical properties of the unzipping 
transition will be modified. In particular, they can yield 
information on the on the three-dimension wandering ex- 
ponent of a single vortex in a disordered sample. 
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B. Unzipping a Luttinger liquid: The clean case 

Consider first an experiment where an attractive two- 
dimensional potential holds vortices confined to it. A 
MFM then pulls a single vortex out of the plane (see 
Fig. [13]). We assume throughout that the density of 
vortices in the three dimensional bulk is so small that we 
can neglect interactions between the vortex that is pulled 
out of the sample and vortices in the three dimensional 
bulk. In this subsection only the clean case (no point 
disorder) will be studied. 

We assume the MFM exerts a force f — fx. As in the 
unzipping experiments discussed above we expect that 
for large forces f > fc the vortex will be completely 
pulled out of the two dimensional slab. Similar to the 
case of the unzipping of a single vortex we write the free 
energy of the vortex as a sum of two contributions. The 
first, Tu{Tm), arises from the part of the vortex that is 
outside the two dimensional slab. The second !Fb{Tm) is 
the change in the free-energy of the vortices that remain 
inside the two dimension slab. As before is the length 
along the r direction which is unbound from the two- 
dimensional slab. The free-energy of the unzipped part 
is clearly identical to that calculated in Eq.[5]or explicitly 

^n(T„0 = -/V,„/27 . (81) 

The calculation of the free-energy, J^6(r,„), is some- 
what more involved. Clearly there is a linear contribu- 
tions due to the length removed from the attractive 
potential of the slab. However, in addition there is an 
extra contribution from the energy of the dislocation, 
^d{Tm), (see Fig. [T21) created in the two dimensional 
vortex array. This contribution to the free-energy, as we 
show below, is non-linear and controlled by the Luttinger 
liquid parameter g. This non-linearity results, near the 
unzipping transition, in a sensitivity of the critical prop- 
erties to the value of g. 

We leave the details of the calculation of the dislocation 
energy to Appendix[X]and present here only the key steps 
of derivation. 

In order to satisfy boundary conditions near the inter- 
face one can use the method of images (see Fig. 
The free energy of this dislocation pair can be calcu- 
lated by standard methods (see details in Appendix [X)) . 
In particular, at large t„i it behaves logarithmically (see 
e.g. Ref. I^): 

= ^ln(T„/ao), (82) 

where oq is the short range cutoff of the order of the 
distance between fiux lines. We note that the free energy 
of the dislocation near the interface ([82]) is one half of 
the free energy of a dislocation pair. 

With the energy of the dislocation in hand we can now 
analyze the properties of the unzipped length near the 
transition using the methods used for analyzing the single 
vortex unzipping experiments. The contributions to the 



free energy are from the unzipped part of the vortex and 
the energy of the dislocation. Collecting all the relevant 
terms, near the transition the free energy is given by 

T 

T{Tm) = J^u(T.m) + ^b{Tm) = eT,„ -|- — ln(T„/ao) . (83) 

The probability of finding a certain value of Tm is then 
given by 

where C is the normalization constant. At the transition 
e = the distribution becomes a pure power law in t,„. 
Therefore, the average value of t„j is very sensitive to 
the value of 5. In particular, for g > 1/4, (i.e. for weakly 
interacting flux lines) the behavior of {Tm) near the tran- 
sition is identical to that of a single vortex in the absence 
of interactions with other vortices 

(t„0 ~ - . (85) 
e 

In contrast, for 1/8 < g < 1/4 (stronger interactions) 
there is a continuously varying exponent governing the 
transition 

(-™) - ^Jrji^ ■ (86) 

And finally, for g < 1/8 (strongly interacting fiux lines) 
we find that (Tm) does not diverge near the transition. 
Note that even though in this regime the mean displace- 
ment remains constant at the transition the higher mo- 
ments of Tm diverge and are thus sensitive to e. The 
reason for this is at the transition the distribution of Tm 
is a power law. 

C. Unzipping from a twin plane with point 
disorder 

We now consider the problem of unzipping a vortex 
from a plane with many vortices in the presence of dis- 
order. In the spirit of the treatments presented in this 
paper, one needs to calculate the free-energy of the un- 
zipped part of the vortex !Fu{Tm), the free-energy of the 
bound part of the vortex !Fb{Tm) and the fluctuations 
in both quantities averaged over realizations of disorder. 
This can be done perturbatively near g = 1. We again 
relegate details of the derivation of the dislocation energy 
to Appendix |B] One conclusion from our calculations is 
that the mean free energy of the dislocation near the 
boundary is not affected by the disorder and is given by 
Eq. ([82]). Another important conclusion is that the fiuc- 
tuations of the free energy also depend logarithmically 



5Tl(Tm) = ln(r„/ao) (87) 
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for g > 1 and 
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In (T,„/ao) (88) 



for g < 1. We note that in the case of many flux hnes 




FIG. 14: Possible experimental setup for studying unzipping 
from Luttinger Liquid in the presence of disorder. 

there is a weak logarithmic dependence of free energy 
fluctuations on as opposed to strong power law depen- 
dence in the case of a single flux line (compare Eq. (|14p V 
This somewhat surprising result is a consequence of the 
screening of strong power-law fluctuations by other flux 
lines. We note that if the pinning of flux lines by disor- 
der is extremely strong so that tearing a single flux line 
does not affect positions of other lines in the duration of 
experiment, we are back to the single flux line physics 
and SJ^J oc t,„. 

To complete the analysis, we need to consider the free- 
energy contribution from the unzipped part. Of partic- 
ular of importance are the free-energy fluctuations due 
to the disorder in the bulk of the sample. As discussed 
in Sec. |TT1 in a three dimensional sample these grow as 
SJ^u oc to'^'^^^ with a;(3) ~ 0.22. This contribution grows 
much quicker than the contribution from the fluctuations 
in the free-energy of the dislocation. Therefore following 
the ideas of Sec. [ill the total free-energy is given by 



Tim) = a(/, - /)t„ - 6t;:;(3) 



(89) 



where a and b are positive constants. Minimizing Eq. 
(|89)) gives for the critical properties in this case 



(90) 



Thus screening disorder fluctuations in the plain by other 
flux lines effectively enhances the role of disorder in the 
bulk. As the result these unzipping experiments can serve 
as a probe of the three-dimensional anomalous wandering 
exponent. 
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APPENDIX A: CALCULATION OF THE 
DISLOCATION ENERGY NEAR THE 
INTERFACE. 

To calculate the energy of the dislocation created by 
unzipping in Fig. 1131 standard methods can be used with 
a slight modification. The need for a modification is due 
to the requirement that vortex lines must exit normal 
to the interface of the superconductor at the top of the 
sample which ensures that the supercurrents are confined 
to the sample. One then has to sum over phonon field 
configuration, Ud{x,T), which satisfy 



dud{x, t) 







(Al) 



T=0 



Here, we have chosen the upper boundary of the two- 
dimensional slab to be located at r = 0. The calculation 
can then be done using the method of images, depicted 
in Fig. [13 by writing the constrained field Ud{x,T) in 
terms of an unconstrained field u{x,t) 



Ud{x, t) = u{x, t) -I- u{x, -t) , 



(A2) 



and summing over all configurations of the unconstrained 
field. In terms of the constrainted field the action is 
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dT I dx[{dxUd{x,T)f + {drUd{x,T)f 



(A3) 



which rewritten in terms of the unconstrained field be- 
comes 



S = ^ 



2g 



dr i dx\{dxu{x,T)y + {dru{x,T))'^ 

! J —OC 

+ dxu{x, T)dxu{x, — r) -I- dru(x, T)dru{x, — r)] (.A4) 
In terms of the partial Fourier transform 
1 



u{x^ t) 



(27r)i/2 



duje"^^u{x,uj) 



(A5) 



the action becomes 



2^/ 



dx [{dxu'{x, uj)Y -t- uj^{u'{x, uj)f] 

) 

(A6) 
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where u'(a;,w) is the real part of u(x,uj). The spatial 
Fourier transform 

oj) = J dqe^^^u'iq, cu) , (A7) 



then finally gives 



OO /'OO 



OO 'J — OO 



S=-l dq I dLo[iq^+LU^)\u'iq,L,)\^] . (AS) 




X 



FIG. 15: An illustration of configurations of the phonon field 
considered by the method of images. 

Eq. IA8I allows a straightforward calculation of the en- 
ergy of a dislocation near the upper boundary of the slab 
(see Fig. [15]) using standard tools. We use the correlation 
function^ » 



G(r„ 



.)] 



(A9) 



where the boson phase angle (f) is conjugate to du/dx. 
Here m is, as before, the length of the unzipped segment 
and the angular brackets denote an average with respect 
to the action Eq. IA8I In terms of this correlation function 
the free-energy of the dislocation is given by 



T 

^d{Tm) = - — In G(t,„) , 



(AlO) 



where T is the temperature and we have set, as before, 
the Boltzmann constant to he ks = 1. Note that we 
are interested in the correlation function of the uncon- 
strained field u{x,t). The factor of one-half originates 
from the fact that the we are interested in the energy up 



to r = and not between — m and m. Integrating out 
the field </> in the standard way it is easy to show that 
this expression can be rewritten as 



G(r) 



exp 



9 J-r„ 

It is then straightforward to find that for large m 

T 



(All) 



Td= -r ln(T„/ao) 



(A12) 



where ao is a microscopic cutoff. The result turns out 
to be identical to that which would have been obtained 
without the constraint on the fields. 



APPENDIX B: CALCULATION OF THE 
DISLOCATION ENERGY IN THE PRESENCE 
OF DISORDER 

As in the clean sample we consider the correlation func- 
tion 



G(T„i) = ( exp - 



9 J-T, 



dxU{0,T)dT 



(Bl) 



Note that in the clean case the method of images guar- 
anteed that the flux lines exit normal to the interface. 
However, with disorder this will be the case only if disor- 
der in the top image plane is correlated with the disorder 
in the bottom plane. However, we expect that at large 
Tm these correlations will not be important and we will 
ignore them. 

To obtain the free-energy with disorder and sample-to- 
sample fl uctuatio ns we consider the replicated correlation 
function G(Tm)". This allows the extraction of the dis- 
order averaged free-energy from 



d 



^d{Tm) 

and the fluctuations in it from 

92 



r^G(r™)"|„=o 



(B2) 



5TI{t^) = r^^G(T™)»|„=o - T' I ^G(t™)"|„=o 



d 



Using standard methods it is straightforward to obtain 
1 



G(r„ 



exp 



2g2 



(B3) 
iain 

dlg{l - Jo{Tme~^/ao)) 



dlK{l){l - Joir^e Voo)) 



, (B4) 



where ao is a microscopic cutoff. The Bessel function 
Jo (a;) appearing in Eq. IB4I and below is non- universal 
and depends on the details of the cutoff procedure. Com- 
bining this with the solutions of the flow equations yields 



Pd{Tm) 



T 



ln(Tm/ao) . 



(B5) 
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Here we have included a factor of one-half since vortices 
exist in the lower half plane. The disordered free energy 
of the dislocation is the same as that of a dislocation in 
the clean sample. Similarly one can obtain for the free- 
energy fluctuations of the dislocation 



for g < 1. 



T 



KlOO 



ln(T„/ao) 



(B6) 



for g > 1 and 



SJ^jiTm) - T'^^—^ln^Tra/ao) (B7) 



As discussed in Sec. IIII CI a crucial input from the 
above calculation is that the free-energy fluctuations 
grow at most logarithmically in Tm- 
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